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Abstract: 

We present the one-loop anomalous dimension matrices for the chiral-even twist-3 (nonsin- 
glet) conformal operators, which govern the scale-dependence of the vector meson light-cone 
wave functions through the conformal expansion. It is clarified that the constraints from 
the charge-conjugation invariance and the chirality conservation allow only one independent 
anomalous dimension matrix for each conformal spin. 
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In a recent paper ffl], we presented a systematic exploration on the twist-3 light-cone 
wave functions (distribution amplitudes) of the vector mesons. These wave functions are 
relevant for understanding preasymptotic corrections to various hard exclusive processes 
producing vector mesons in the final state |2], [3] . According to that study, two chiral-even 
wave functions g± (u), g± (u), and two chiral-odd ones hi? (it), hi (u) constitute a complete 
set of the twist-3 quark- ant iquark wave functions (u denotes the light-cone momentum 
fraction carried by quark). The renormalization of these wave functions has also been studied 
based on the (approximate) conformal invariance in massless QCD Q: The wave function 
can be decomposed into "partial waves" of definite conformal spin, and the coefficients in 
this expansion are given as vacuum-to-meson matrix elements of the conformal operators. 
Conformal symmetry allows renormalization mixing only among the conformal operators 
with the same conformal spin up to 0(a 2 ) correction. It has also been shown in |IJ that the 
anomalous dimensions of those conformal operators can be expressed in terms of those for the 
nucleon's structure functions (parton distributions) for inclusive processes. In particular, the 
"/i 2 -evolution" of the chiral-odd wave functions hi? (u, /x 2 ) and hi (u, /i 2 ) has been solved 
in terms of the anomalous dimension matrices for the nucleon's chiral-odd twist-3 quark 
distribution functions e(x, /i 2 ) and /i^(a;,/i 2 ) || |6], [7], §. On the other hand, the anomalous 
dimension matrices for the chiral-even wave functions have been only partially given in ]l| 
based on the known results for the nucleon's another chiral-even distribution, grp(x,n 2 ) |J 
0. 

The purpose of this letter is to complete the renormalization of the twist-3 conformal 
operators relevant to chiral-even wave functions g±(u,/i 2 ) and g± \u, /i 2 ). We give a con- 
venient choice of the independent operator basis, and clarify the symmetry constraints on 
their renormalization mixing. We shall complete the relevant anomalous dimension matrices 
for all conformal spins. 

Let us first introduce a twist-3 conformal operator basis relevant for the renormalization 
of the nonsinglet chiral-even wave functions, g±\u, fi 2 ) and g±\u, fi 2 ), and give precise 
formulation of our problem referring to the results of M. The QCD equations of motion allow 
us to reexpress the quark-antiquark wave functions g±\u, jj, 2 ) and g± (u, /i 2 ) in terms of the 
twist-3 quark-antiquark-gluon wave functions. This can be performed explicitly order by 
order in conformal expansion. For the conformal partial wave of conformal spin j = n + 3/2 
(n = 2,3, •••), two sets of twist-3 quark-antiquark-gluon conformal operators, {-R^ fc } and 
come into play: 

Kn,k = ®k,n-k-2 ^ ®n-fc-2,fc T ®fc,n-fc-2 + ®n-k-2,k (^ = 0, 1, . . . , 72 — 2) , (1) 

e^ n _ jfc _ 2 (0) = ip(0)d n - k - 2 figG LX (0)A x d k ip(0) + (total derivatives), (2) 

fc,n-fe-2(°) = ip(0)d n - k - 2 figG LX (0)A x d k i^(0) + (total derivatives), (3) 

where we introduced an auxiliary light-like vector A M with A + = A -1 = 0, which guarantees 
the symmetrization of the n Lorentz indices and the subtraction of the trace terms, and 
d ee iA-D with = d^ — igA^ the covariant derivative. For the purpose of renormalization, 
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it suffices to consider the flavor- diagonal operators as given above. The dual gluon field 
strength is defined as G^ u = ^e^\ a G Xa where e^Ao- is the totally anti-symmetric tensor 
with eoi23 = 1- (We follow the convention of [[I in this paper). In (||), 'total derivatives' 
stands for the terms consisting of 

(A • d) n ~ r [^d r - j -%gG ±x A x d j ip] , ( 2 < r < n - 1, < j < r - 2), 

and similarly for those in (|3|). These total derivative terms are organized so that the opera- 
tors (0) and (Q) form irreducible representation of conformal spin j of the so-called collinear 
conformal groupj|]. The complete expression for {O^'n-k-2} * s given in terms of the Ap- 
pell polynomials, but it is irrelevant here. The suffix k of ([T])-(|3]) labels many independent 
conformal operators having the same conformal spin, illustrating that three-particle repre- 
sentations of conformal group are degenerate; the number of independent operators increases 
with the spin. Conformal symmetry does not allow renormalization mixing in the leading 
logarithmic order between the operators with different n; it does allow, however, mixing 
between the operators with different k for the same n. Here we note that there exists a con- 
straint from symmetry: {R^k} an d {Rnk} have a definite charge-conjugation parity (— l) n 
and (— l) n+1 , respectively, and do not mix with each other under renormalization. If we 
define the anomalous dimension matrix for {R^k\ as T^, the renormalization group (RG) 
equation for {R^ k ] is given by 



a 
2n 



n-2 



(4) 



where fij- = ji-^- + /3(g)-j^ with (3(g) the /5-function. This equation is solved to give the 
scale-dependence of R^ k as 

Kk(0;Q 2 ) = £ \(^WX n/b ] <<(0;^ 2 ), 



(5) 



k,l 



where b 



yiV c — \Nf. The mixing matrix is subject to our computation. The nucleon's 
forward matrix elements of R^ k with even n is associated with the n-th moment of the trans- 
verse spin structure function g2(x,Q 2 ) in the deep inelastic scattering. In this connection 
the renormalization of R^ k has been studied by several different approaches ||, [UJ for even 
n, while the result for odd n has not been considered seriously. On the other hand, there 
have been no places in inclusive processes where R~ k appear, and their renormalization has 
never been discussed before. 



In order to obtain the anomalous dimension matrix for the conformal operator basis 
{Rnk\i it i s convenient to work with their forward matrix elements: We imbed {R^k} 
the three-point function as 



d xd yd ze n 



ipx—iqy+i{q—p)z 



<0|Ti?± m(xMy)A c (z)\0). 



(6) 
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Then, the total derivative terms in (0) and (|3[) drop out and the calculation becomes greatly 
simplified. For later convenience, we introduce two sets of operators, which replace k in 



where k 



0,1,. 



Note that and W%£~~ agree with R^ k and R n k, respectively, if one ignores the 

total derivatives from the latter. For the renormalization, we follow the method of H [7], |T0f : 
We calculate the one-loop corrections for ([]). To renormalize the composite operators 
involving massless fields consistently, we keep the quark and gluon external lines off-shell. 
We then find that the following twist-3 operators also participate: {W^} mix withP] 



W. 



2, and 



ipd 



?t, k 



■n-k-2. 



u (+)k 



° \-)n-k-2i 



4L(gG° x A x ±z l5 gG° x A x ) dfy. 



(7) 
(8) 

(9) 



W. 



(a)o- 



n,F 



w. 



(a)a 



n+1 



n 



k=i 



n,E 



2(n+ 1) 



while {W^k*} niix with 



(10) 
(11) 



ir, 



(v)a 



n,F 



n+1 ^ 



n 



k=i 



n—k 



(12) 



IF, 



(y)o 



n 



n,E 



2( n + l)* e ^ A 



ipiiPY^d n - 1 ip + ipYfid^Hl/)^ 



(13) 



where we introduced another light-like vector £ M with t~ — i — and i ■ A = 1. In the 
following, we assume and A M have mass dimensions +1 and —1, respectively. The for- 
ward matrix elements of the quark bilinear operators W^fp and W^p coincide with those 

of the twist-3 quark-antiquark conformal operators of spin j = n + 3/2 for g^\u, fi 2 ) and 
g^L (u, fi 2 ). It is instructive to note that they originate from the light-cone nonlocal operators 
^(0)757 <T [0, A]<0(A) and ^(0)7 a [0, A]^(A), where [0,A] = exp {-ig dtA^(tA)} is the 
gauge link operator: Taylor expansion of these nonlocal operators at small quark-antiquark 
separations yields ^("f^^iD^ 1 ■ ■ -iD^ip, and the subtraction of the twist-2 components, 

2 The operators Wj^p, 2 (n+i) ^nfe* 7 ' an( l ^n"lT coincide with a set of operators considered in jlo| for 
the massless quark limit. 
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which are totally symmetric and traceless among {a, fii, ■ ■ ■ , fi n }, leads to W^p a corre- 
sponding to antisymmetric pairs of indices a and W^g are the so-called equation- 
of-motion (EOM) operators, which contribute as nonzero operators to the off-shell Green's 
function. As it is well known, the above twist-3 operators are not independent, but obey 
the following relation: 

-i n-2 

wiy h = j> + i)w^ + (i4) 

l \ n + L ) k=0 

Therefore, we can conveniently choose (0) and ( p]) ((§) and (|T3"D) as an independent basis 
for renormalization. A physical matrix element of vanishes, and does not affect the 

final results of the scale-dependence. However, it has to be kept in the intermediate step of 
the calculation to work out the renormalization [|, 0, [U| . 

Here we comment on the places where these operators appear in inclusive processes in- 
volving nucleons. The nucleon matrix element (PS\W^ p (Q 2 )\PS) gives the twist-3 part 
of the n-th. moment j\ dxx n gT{x, Q 2 ) with gx(x : Q 2 ) the nucleon's transverse parton dis- 
tribution. (| -PS 1 ) is the nucleon state with momentum P M and spin vector S^). For even n, 
this moment is associated with Jq 1 dxx n g2(x, Q 2 ) measurable in the deep inelastic scattering, 
while both even and odd moments are relevant in the polarized Drell-Yan process (see e.g. 
fnj). In order to identify the role of W^p in inclusive processes, we recall that it is gen- 



erated from the nonlocal light-cone operator , 0(O)7 cr [O, A]^(A) and consider the following 
decomposition of the nucleon matrix element, 

r d\ 

j — e lXx (PS\m^, XA]^(XA)\PS) = 2 [f^x)^ + S a J T (x) + M 2 A fl f 4 (x 

(15) 

where P» and S» are decomposed as P" = P + ^-A^ and S" = (S ■ A)^ + (S ■ £)A» + 
with the nucleon mass M and P 2 = —S 2 = M 2 . fi(x) and fi(x) are the spin- independent 
parton distribution functions of twist-2 and twist-4, respectively, as discussed in [fLl}| . The 
new function fr(x) is of twist-3 112[ . By comparing both sides of (|15]), it is easy to see that 
j\dxx n fT(x,Q 2 ) is associated with (PS\wj l p(Q 2 )\PS}. However, the existence of fr(x) 

violates time reversal invariance, and therefore (PSIW^PS) = (PSIW^PS) = 0. 
(The nonforward matrix elements like (0\R^ k \P) relevant to meson wave functions are not 
constrained to vanish by time reversal invariance.) 

In order to work out the renormalization of W^^ a , we clarify the symmetry constraints 
on the renormalization mixing. We find: 

(i) Two sets {U? + j k } and {U7_s k } in @ do not mix with each other. To see this, we 
recall o =_L and introduce right- and left-handed circular polarization vectors for the 
gluon as R = (0, 0)/v2 and L = (0, 1, — i, 0)/y/2, respectively. Then we find 
U(=^ = -2i^) L d n - k - 2 ^G Rv 'A v d k ip L and = -2%i) R d n - k - 2 l^G Rh 'A v d k ip R where 

i>R,L = (l±75)^/2. We also obtain similar expressions with R <-» L for Uh^ k and Up^ k . 
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These results explain no mixing between {Uf + ^ k } and {U?\ k }, since the perturbative 
quark- gluon coupling preserves chirality of each quark line. 

(ii) The anomalous dimension matrices for {U^ k } and {U? \ k } are identical: From (i), 
the RG equation for {Uf^ k } can be written with its anomalous dimension matrix T n 

as 

= -lfgVj^<v <«>> 

where the terms due to the EOM operators are ignored assuming that we take a 
certain physical matrix element. By applying the charge conjugation C on both sides, 
one obtains 

d a s n ~ 2 

A t ^^(-)n-fc-2 = ( r ™)fe,/ ^f-)ti-l-2! (17) 

where we used CU^ k C~ 1 = (—l) n U^ n _ k _ 2 . Equations (|T6|) and fllTD explain the 
statement. 

From fll6|) and (]T7D , one obtains the RG equation for as 

;7 n— 2 



By comparing this result with forward matrix element of (§]), we can identify the anomalous 
dimension matrix as 

(pn ) fc l = (Fn)n-k-2,n-l-2 i (-^n ) fe ; = C^n)^ • (19) 

We also note that we can use whichever convenient set of the operators {W^ k ^ a }, {U^ k }, in 
order to obtain T n (see (p~^) - (p~8|) ) . These features of the anomalous dimensions are in parallel 
with a familiar example of the twist- 2 distributions fx and g\ defined from the nonlocal 
operators ^(0)7 CT [0, A]-0(A) and ^>(0)757 CT [0, A]^(A): Their n-th moment is associated with 
the local operators ipd n ^Lip = tjjRd n /^ipR + ^jjLd n /^^jjL and -0<i n ^75^ = ^^"A^h — ^lgPA^l, 
respectively, and similar argument as above explains that the nonsinglet parts of f\ and g± 
have the same anomalous dimensions. 



Explicit calculation of the mixing matrix r„ is yet desirable for general (even and odd) n. 
As noted above, one has to keep the EOM operators as nonzero operators during the course 
of the renormalization. We choose {W^ k ^ a , W^ a ^ a } as a basis of the renormalization. (For 
completeness, we parallely discuss two sets of the operators.) We imbed these operators 
into the three-point function (^). The relevant one-loop diagrams are the same as those 
considered in [|10j. We adopt the Feynman gauge, and obtain the mixing matrix T n in the 
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MS scheme of the dimensional regularization. One technical difficulty in the computation 
is the complicated mixing with many gauge noninvariant EOM operators. To avoid this, 
we introduce a vector with the condition Q ■ A = 0, and contract the external gluon line 
of ^ 

written as 



with this vector 0. Under this condition, the vertices for W^i^'* and W^ a F >a can be 



r(a,V)(T 



W { :f a = ±g(q - p) (p fc g- fe - 2 757W =F fr-^W^) t c , 



n 
+ 



(20) 
(21) 



where t c is the color matrix normalized as tr(t c t c ) = \5 CC \ and we introduced the notation 
p — p ■ A, etc. It is now straightforward to obtain the mixing matrix T n . The result reads 

(/ + !)(/ + 2) 



(T 



nJkl 



-a 



G 



(k + l)(k + 2)(k-l) 

k+i+i n-iQ n + k — I 



(T 



+(2C F - C G 
-1 



-1 



nJkk 



a 



G 



1 



-iCfc n(k - 1) 
1 1 



+ (-1 



n — k — 1 n — k n — k + 1 k + 2 

-1 2(-ir- fc 



+ (2C F -C G 



k + 2 

+ 'S'n-fe-l + Sk+1 



(k > I) (22) 



n 



(r 



(n — k — l)(n — A;)(n — A; + 1) 
+C*i? (25*fc + i + 2S' n _/ c _i — 3) , 
(n - l)(n - I + 1) 



k + 2 



(23) 



-a 



G 



(n - k)(n - k + - k) 



+ (2c F - a 



G 



-i 



k+l+ 



i n-iQ+i n + l-k 



2{-iy 



n—k 



-Pi- 



(n — k — l)(n — k)(n — k + 1) 



n-Pk+i n{l-k) 
(I > k) (24) 



where n C k = n\/[k\{n - k)% S n = ELi and C F = {N 2 C - 1)/2N C , C G = N c for 

N c color. This result combined with (|19|) gives T^. The result for T~ is wholly new. For 
T^, the result is now completely applicable also for the cases with odd n which has not 
received much attention.^ For example, if one replaces the term (— 1) / (n — I) proportional 
to 2C F — C G in Eq. (16) of [|1(J by (— l)"~ l+1 /(n — /) in their convention, their result in 
Eqs. (15)-(17) becomes consistent with the above (|22"|) -(f2"4"l) for all n. 

We make contact with the work of Balitsky and Braun|13[, who also considered the 
renormalization of similar quark-antiquark-gluon operators. They introduced the nonlocal 



3 T 2l T+, and IV give Eqs. (4.66), (4.67) and (4.68) of ||. 
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light-cone operators 



#(aA)[aA,/?A] (G W (J3A) ± G av {ftA)i l5 ) A v A[ftA, T A]^( T A), (25) 



and computed the one-loop evolution kernels in the coordinate space representation. They 
found by explicit calculation that S^(a,ft,"f) and S~(a,ft,j) are renormalized without 
mixing with each other. The reason for this phenomenon is made clear by examining the 
R, L components of similarly to (i) above. A direct relation between our results 



a 



and those of |TB| can be provided by expanding the operators ( P5|) in powers of A u and 
taking the forward matrix element: Expansion of S^{a, ft, 7) generates a series of the local 
operators with increasing spin. In particular, the n-th term generates A u £ x e au \^U^ ± - )k with 
k = 0, 1, ... , n— 2, up to the twist-4 corrections. In this connection, we mention the following 
point: The kernels for S£(a,ft,j) and S~(a, ft, 7), given by Eq.(6.2) of have rather 
different forms, although we have shown that the anomalous dimension matrices for the 
corresponding two sets of local operators, {U7 + ^ k } and {U^ k }, are identical with each other 
(see ( |I6D and flT7D). In fact, the equivalence of these two kernels can be easily proved by 
applying the parity transformation combined with the time-reversal transformation^ to the 



corresponding evolution equations, and therefore the results of |L3[ are consistent with ours. 

We finally recall in our notation that the anomalous dimension matrix in (|22"|)-(|2"4l 
satisfies the relation in the large N c limit, i.e. Cp — > N c /2 in (^2|) -(^) fl4]| : 



n-2 



E(*+i)(r 



(I + l)7n, 



(26) 



fc=0 



where 7„ is the lowest eigenvalue of the anomalous dimension matrix T n in this limit, and 
is given by 



ln = 2N e L(n + l) +lE -\ + — l — 



(27) 



Here ip(n + 1) = J2k=i 1/^ — 7s is the digamma function and 7^ is the Euler constant. This 
leads to a simple evolution equation in the large N c limit (see (||), (0), and (|T9|)): 



a(Q 2 ) 



a s (/i 2 



ao^ 2 ), 



(28) 



where C„(/i 2 ) = £ 



fc=0 



n- 



■ fc - l)< fe (0; /i 2 ), E^=o 2 (^ + 1)^(0; A* 



2 \ and W nF 



(a,v)a i 



(/^ 2 ). Similar 

relations ( p6|) and (p8|) hold also for the limit of large conformal spin with j n — > 7„ + (4Cj? — 
2iV c )(lnn + 7^ — 3/4) to 0(\nn/n) accuracy. As a result of (EH), the vector meson wave 



functions <7j°' '(«, /i 2 ), g±'(u, /i 2 ) as well as the nucleon parton distribution g T (x,jj 2 ) obey 
simple DGLAP-type evolution in the two limits, iV c — > 00 and n — > 00 fL], |14 



This transformation is more convenient for this purpose than the charge-conjugation transformation 
used in ( |l6| ) and (|l7|), since it connects S£{a,(3, 7) with S~(a,P,j) keeping the condition a > 7, which is 
assumed in Eq.(6.2) of [O. 
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To summarize, we have completed the renormalization of chiral-even twist-three non- 
singlet wave functions for the vector mesons. The result is presented in the form of the 
anomalous dimension matrices for the corresponding conformal operators {R^k) = 
0, 1, • • • , n— 2) for all conformal spins. Although {R^k} an d {Rnk} are the independent oper- 
ators, their anomalous dimension matrices coincide completely due to the charge-conjugation 
invariance and the chirality conservation. 
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